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(54) Method and apparatus for automatic setting of sample gate in pulsed doppler ultrasound 
imaging 



(57) A method and an apparatus for automatically 
initializing and adjusting the Doppler sample gate posi- 
tion and size settings based on actual vessel image da- 
ta. A vessel segment search method employs an object 
search technique based solely on geometric and mor- 
phological information in a binarized vessel image ob- 
tained from either B-mode or color flow image data. The 



morphologically best or nearest vessel segment within 
a target search region in the two-dimensional image is 
found. The sample gate is placed at or near the center 
of the targeted vessel segment. The sample gate size 
is adjusted in relation to the vessel size. Then the best 
available steering angle that minimizes the Doppler an- 
gle is selected. 
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Description 

[0001 ] This invention generally relates to the imaging 
of moving ultrasound scatterers. In particular, the inven- 
tion relates to methods for positioning the gate or sam- 
ple volume (hereinafter "sample gate") in medical diag- 
nostic ultrasound imaging. 

[0002] Premium medical diagnostic ultrasound imag- 
ing systems require a comprehensive set of imaging 
modes. These are the major imaging modes used in 
clinical diagnosis and include timeline Doppler, color 
flow Doppler, B mode and M mode. In the B mode, such 
ultrasound imaging systems create two-dimensional im- 
ages of tissue in which the brightness of a pixel is based 
on the intensity of the echo return. Alternatively, in a 
color flow imaging mode, the movement of fluid (e.g., 
blood) or tissue can be imaged. Measurement of blood 
flow in the heart and vessels using the Doppler effect is 
well known. The phase shift of backscattered ultrasound 
waves may be used to measure the velocity of the back- 
scatterers from tissue or blood. The Doppler shift may 
be displayed using different colors to represent speed 
and direction of flow. In the spectral Doppler imaging 
mode, the power spectrum of these Doppler frequency 
shifts are computed for visual display as velocity-time 
waveforms. 

[0003] One of the primary advantages of Doppler ul- 
trasound is that it can provide noninvasive and quanti- 
tative measurements of blood flow in vessels. Given the 
angle 6 between the insonifying beam and the flow axis 
(hereinafter referred to as the "Doppler angle"), the 
magnitude of the velocity vector can be determined by 
the standard Doppler equation: 

v=c/ d /(2/ 0 cos9) (1) 

where c is the speed of sound in blood, / 0 is the transmit 
frequency and / d is the motion-induced Doppler fre- 
quency shift in the backscattered ultrasound signal. 
[0004] In conventional ultrasound spectral Doppler 
imaging, the operator is required to manually position 
the sample gate to the measurement location in a two- 
dimensional image with or without color flow data. The 
operator also needs to manually adjust the sample gate 
size relative to the diameter of the vessel to be studied. 
From the acoustic data acquired over many transmit fir- 
ings, Doppler frequency spectral data is obtained via 
standard Fast Fourier Transform (FFT) spectral analy- 
sis. 

[0005] For a given measured Doppler frequency shift 
/ d , the flow velocity (speed) v is calculated using Eq. 
(1 ). The ideal Doppler angle is zero, i.e., when the beam 
is aligned in the direction of blood flow. Unfortunately, 
Doppler angles that can be formed in practice tend to 
be larger, and as they approach 90 degrees, a small er- 
ror in the angle estimate can lead to a large error in v. 



2 

For this reason, it is generally recommended that the 
Doppler beam be steered to form Doppler angles of no 
greater than about 60 degrees for reliable velocity 
measurements. Angle steering is another adjustment 
5 the operator needs to make manually in conventional 
Doppler systems. 

[0006] In an attempt to minimize manual Doppler ad- 
justments when the spectral Doppler mode is activated, 
conventional scanners generally provide presets for the 

10 Doppler sample gate position and size, and forthe beam 
steering angle. However, such presets have limited ben- 
efits because the vessel depth, size and orientation rel- 
ative to the probe can vary a great deal from one case 
study to the next. 

15 [0007] U.S. Patent No. 5,365,929 describes the use 
of multiple range gates and multiple Doppler beams to 
scan a region of interest. By comparing some signal 
characteristic, such as total power or maximum velocity, 
of the multiple sample volumes, the scanner automati- 

20 cally selects the best sample gate for full spectral anal- 
ysis and display. It will appear to the user that the scan- 
ner has automatically positioned the sample gate at a 
location where the Doppler signal is optimal in some 
sense. 

25 [0008] European Patent Application No. 0 842 638 A2 
describes a method of tracking vessel walls in the B- 
mode image, and then automatically adjusting the sam- 
ple volume size to ensure the entire vessel diameter is 
covered. While this may be useful for volume flow meas- 

30 urements, the user is still expected to first manually po- 
sition the sample volume and vessel wall markers at the 
correct locations. Also, in Doppler exams that do not in- 
volve volume flow measurements, different clinics may 
follow different practices in terms of sample gate size 

35 relative to the vessel diameter. 

[0009] European Patent Application No. 0 985 380 A1 
describes a method for automatic positioning of the 
Doppler sample gate based on bloodstream or color 
flow information. Among various specific applications, 

40 this method can be used to automatically set the sample 
gate cursor at an optimal position when the sample gate 
is first brought up in the image, or when it is being 
moved. The optimal position may be defined by a color 
flow pixel showing the highest velocity, or the center 

45 point of the largest flow segment, or the center point of 
the next best flow segment etc. 
[001 0] U.S. Patent No. 5,690,1 1 6 describes a method 
of estimating the orientation (slope) of the vessel axis 
based on gray-scale image data, and then computing 

so the Doppler angle. 

[0011] U.S. Patent No. 6,068,598 describes a robust 
method for detecting the vessel walls based on B-mode 
and/or color flow data, estimating the vessel orientation 
based on the best vessel edge data, and then computing 

55 the Doppler angle. 

[0012] U.S. Patent No. 4,937,797 describes a method 
of adjusting the transducer array beam-forming delays 
to automatically steer the Doppler beam to achieve a 
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target Doppler angle such as 60 degrees (or less). This 
method, however, requires the user to first manually ro- 
tate an angle cursor on the B-mode image to define the 
vessel orientation or flow direction. 
[0013] There is a need for an automatic method of in- 
itializing and adjusting the Doppler sample gate position 
and size settings and the beam steering angle setting 
based on actual vessel image data, with the goal of im- 
proving the efficiency of the Doppler study above and 
beyond what can be achieved using presets. 
[0014] The present invention is directed to a method 
and an apparatus which provides automatic initialization 
and adjustment of the Doppler sample gate position and 
size settings, and of the beam steering angle setting 
based on actual vessel image data. This capability is 
presented, in one embodiment, as an aid or "smart fea- 
ture" to the user when either one of the two following 
conditions is met: (1) when the sample gate cursor is 
first activated or brought up in the image (before or upon 
spectral Doppler mode activation); or (2) when the sam- 
ple gate cursor is being moved. Both of these conditions 
are normally monitored automatically by the master con- 
troller of the system. In accordance with the preferred 
embodiment of the invention, if either condition is de- 
tected, the following functions are automatically per- 
formed: (1 ) the "morphologically best" or nearest vessel 
segment within a target search region in the two-dimen- 
sional image is found; (2) the sample gate is placed at 
or near the center of the targeted vessel segment; (3) 
the sample gate size is adjusted in relation to the vessel 
size; and (4) the best available steering angle that min- 
imizes the Doppler angle is selected. All four features 
are intended to provide a good but initial adjustment to 
the sample gate position and size and to the beam steer- 
ing angle, which should facilitate further manual adjust- 
ments by the user. It is not generally possible to predict 
exactly which vessel or vessel portion the user wants to 
study, but by setting up the Doppler sampling parame- 
ters according to the best image data available, this in- 
vention will improve examination speed and/or ease of 
use. 

[0015] The vessel segment search method (step 1 
above) in accordance with the preferred embodiment of 
the invention employs an object search technique based 
solely on geometric and morphological information in a 
binarized vessel image obtained from either B-mode or 
colorf low image data. The vessel extraction method dis- 
closed herein is not based on comparing Doppler signal 
characteristics from multiple sample volumes. The 
method in accordance with the preferred embodiment 
comprises the following steps: (a) reading an image 
frame from memory; (b) defining a search region within 
the image frame or a decimated version of the image 
frame; (c) binarizing the pixel values within the search 
region based on a pixel value threshold that is adaptive 
to local pixel value statistics or based on the presence 
or absence of color flow information in each pixel; (d) 
morphological filtering the binarized pixel values to elim- 



inate structures smaller than the speckle size; (e) count- 
ing the number of continuous objects and keeping track 
of their area; (f) rejecting objects that are outside a pre- 
determined size range; and (g) finding the center point 

5 of the "best" vessel segment. As a result of this method, 
the sample gate position can be automatically set. 
[001 6] Having identified the center of the best vessel 
segment, the vessel diameter and orientation can be 
computed using any automated method. The sample 

10 gate size can then be automatically set to be some frac- 
tion of the vessef diameter. 

[001 7] The invention will now be described in greater 
detail, by way of example, with reference to the draw- 
ings, in which:- 

15 

FIG. 1 is a schematic showing a block diagram of a 
typical ultrasound imaging system which can be 
programmed with software in accordance with the 
preferred embodiment of the present invention. 

20 

FIG. 2 is a schematic depicting an ultrasound image 
of a portion of a blood vessel with a sample gate 
graphic, a Doppler beam cursor and a vessel slope 
cursor superimposed thereon. 

25 

FIG. 3 is a schematic depicting an ultrasound image 
with a search region for sample gate re-selection 
shown as a dashed rectangle in one quadrant. 

30 FIG. 4 is a flowchart showing an algorithm for auto- 
matic detection of vessel segments in accordance 
with the preferred embodiment of the invention. 

FIG. 5 is a flowchart showing an algorithm for bina- 
35 rization of the search region in accordance with the 
preferred embodiment of the invention. 

FIG. 6 is a schematic depicting an ultrasound image 
of a portion of a blood vessel with edge point search 
40 information superimposed thereon. 

FIG. 7 is a schematic depicting segmentation of the 
edge points in accordance a method for automated 
estimation of the Doppler angle. 

45 

FIG. 8 is a flowchart depicting the algorithm for au- 
tomatically adjusting the vessel slope cursor in ac- 
cordance with a method for automated estimation 
of the Doppler angle. 

50 

[0018] An ultrasound imaging system programmed 
with software in accordance with the preferred embod- 
iment of the present invention is generally depicted in 
FIG. 1 . The main data path begins with the analog RF 
55 inputs to the beam-former board 4 from the transducer 
2. The beam-former board 4 is responsible for the trans- 
mit and receive beam-forming. The beam-former's sig- 
nal inputs are the low-level analog RF signals from the 
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transducer elements. The beam-former board 4, which 
comprises a beam-former, a demodulator and filters, 
outputs two summed digital baseband I and Q receive 
beams formed from acquired data samples. These data 
samples are derived from the reflected ultrasound from 
respective focal zones of the transmitted beams. The I 
and Q data is sent to FIR filters which are programmed 
with filter coefficients to pass a band of frequencies cen- 
tered at the fundamental frequency f 0 of the transmit 
waveform or a (sub)harmonic frequency thereof. 
[001 9] The image data output from the filters is sent 
to the midprocessor subsystem, where it is processed 
according to the acquisition mode and output as proc- 
essed vector data. Typically, the midprocessor subsys- 
tem comprises a color flow processor 6, a B-mode proc- 
essor 8 and a spectral Doppler processor 10. Alterna- 
tively, a digital signal processor or array of such proces- 
sors can be programmed to process signals for all three 
modes. 

[0020] The B-mode processor 8 converts the base- 
band I and Q data from the beam-former board 4 into a 
log-compressed version of the signal envelope. The B- 
mode function images the time-varying amplitude of the 
envelope of the signal as a gray scale. The envelope of 
a baseband signal is the magnitude of the vector which 
I and Q represent. The l,Q phase angle is not used in 
the B-mode display. The magnitude of the signal is the 
square root of the sum of the squares of the orthogonal 
components, i.e., (P+CP) V2 .lhe B-mode intensity data 
is output to a B-mode acoustic line memory (not shown) 
in the scan converter 12. 

[0021] The scan converter 12 accepts the processed 
B-mode vector data, interpolates where necessary, and 
converts the data into X-Vformat for video display. The 
scan-converted frames are passed to a video processor 
1 4, which maps the video data to a gray-scale mapping 
for video display. A conventional ultrasound imaging 
system typically employs a variety of gray maps, which 
are simple transfer functions of the raw image data to 
display gray levels. The gray-scale image frames are 
then sent to the display monitor 1 6 for display. 
[0022] The B-mode images displayed by monitor 16 
are produced from an image frame of data in which each 
datum indicates the intensity or brightness of a respec- 
tive pixel in the display. An image frame may, e.g., com- 
prise a 256 x 256 data array in which each intensity 
datum is an 8-bit binary number that indicates pixel 
brightness. Each pixel has an intensity value which is a 
function of the backscatter cross section of a respective 
sample volume in response to interrogating ultrasonic 
pulses and the gray map employed. The displayed im- 
age represents the tissue and/or blood flow in a plane 
through the body being imaged. 
[0023] The color flow processor 6 is used to provide 
a real-time two-dimensional image of blood velocity in 
the imaging plane. The frequency of sound waves re- 
flecting from the inside of blood vessels, heart cavities, 
etc. is shifted in proportion to the velocity of the blood 



cells: positively shifted for cells moving towards the 
transducer and negatively for those moving away. The 
blood velocity is calculated by measuring the phase shift 
from firing to firing at a specific range gate. Instead of 
5 measuring the Doppler spectrum at one range gate in 
the image, mean blood velocity from multiple vector po- 
sitions and multiple range gates along each vector are 
calculated, and a two-dimensional image is made from 
this information. Thecolorflow processor 6 receives the 
10 summed left and right, complex l/Q data from the beam- 
former board 4 and processes it to calculate the mean 
blood velocity, variance (representing blood turbulence) 
and total pre-normalization power for all sample vol- 
umes within an operator-defined region. These three 
15 output values are then combined into two final outputs, 
one primary and one secondary. The primary output will 
be either velocity or power. The secondary output can 
be either variance or power. Which two values will be 
displayed is determined by the operator-selected dis- 
play mode. Both values are sent to a color acoustic line 
memory (not shown) in the scan converter 12. The 
change or shift in backscattered frequency increases 
when blood flows toward the transducer and decreases 
when blood flows away from the transducer. Typically, 
color flow mode displays hundreds of adjacent sample 
volumes simultaneously, all laid over a B-mode image 
and color-coded to represent each sample volume's ve- 
locity. 

[0024] In the color flow mode of the conventional ul- 
trasound imaging system, an ultrasound transducer ar- 
ray is activated to transmit a series of multi-cycle (typi- 
cally 4-8 cycles) tone bursts which are focused at the 
same transmit focal position with the same transmit 
characteristics. These tone bursts are fired at a pulse 
repetition frequency (PRF). The PRF is typically in the 
kilo-hertz range. A series of transmit firings focused at 
the same transmit focal position are referred to as a 
"packet". Each transmit beam propagates through the 
object being scanned and is reflected by ultrasound 
scatterers such as blood cells. The return signals are 
detected by the elements of the transducer array and 
then formed into a receive beam by a beam-former. 
[0025] For example, the traditional color firing se- 
quence is a series of firings (e.g., tone bursts) along the 
same position, which firings produce the respective re- 
ceive signals: 

F 1 F 2 F 3 F 4 - F M 

where F/ is the receive signal for the Mh firing and M is 
the number of firings in a packet. These receive signals 
are loaded into a comer turner memory, and a high pass 
filter (wall filter) is applied to each down range position 
across firings, i.e., in "slow time". In the simplest case 
of a (1 , -1) wall filter, each range point is filtered to pro- 
duce the respective difference signals: 
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(F r F 2 )(F 2 -F 3 )(F 3 -F 4 ) ,.(F MV F M ) 

and these differences are input to a color flow velocity 
estimator. Typically, the corner turner memory, wall filter 
and parameter (e.g., velocity) estimators are incorporat- 
ed into the color flow processor 6. 
[0026] The color and B-mode acoustic line memories 
in scan converter 12 respectively accept processed dig- 
ital data from the color flow and B-mode processors. 
These components of the scan converter also perform 
the coordinate transformation of the color flow and B- 
mode data from polar coordinate (fl-9) sector format or 
Cartesian coordinate linear format to appropriately 
scaled Cartesian coordinate display pixel data, which is 
stored in an X-Y display memory (not shown) in the scan 
converter. Color flow images are produced by superim- 
posing a color image of the velocity of moving material, 
such as blood, over a black and white anatomical B- 
mode image. 

[0027] If the image to be displayed is a combination 
of one B-mode frame and one color flow frame, then 
both frames are passed to the video processor 1 4, which 
maps the B-mode data to a gray map and maps the color 
flow data to a color map for video display. In the final 
displayed image, the color pixel data is superimposed 
on the gray-scale pixel data. Successive frames of color 
flow and/or B-mode data are stored in a cine memory 
24 on a first-in, first-out basis. 

[0028] Storage can be continuous or as a result of an 
external trigger event. The cine memory 24 is like a cir- 
cular image buffer that runs in the background, captur- 
ing image data that is displayed in real time to the user. 
When the user freezes the system (by operation of an 
appropriate device on the operator interface 22), the us- 
er has the capability to view image data previously cap- 
tured in cine memory. 

[0029] In spectral Doppler imaging, the l/Q compo- 
nents are integrated (summed) over a specific time in- 
terval and then sampled by the spectral Doppler proc- 
essor 10. The summing interval and the transmit burst 
length together define the length of the sample volume 
(i.e., sample gate) as specified by the user. A "sum and 
dump" operation effectively yields the Doppler signal 
backscattered from the sample volume. The Doppler 
signal is passed through a wall filter which rejects any 
clutter in the signal corresponding to stationary or very 
slow-moving tissue. The filtered output is then fed into 
a spectrum analyzer, which typically takes Fast Fourier 
Transforms (FFTs) over a moving time window of 32 to 
1 28 samples. Each FFT power spectrum is compressed 
and then output by the spectral Doppler processor 1 0 to 
a graphics/timeline display memory 1 8. The video proc- 
essor 14 maps the compressed spectral Doppler data 
to a gray scale for display on the monitor 1 6 as a single 
spectral line at a particular time point in the Doppler ve- 
locity (frequency) versus time spectrogram. 
[0030] System control is centered in a host computer 



(i.e., master controller) 20 : which accepts operator in- 
puts through an operator interface 22 (e.g., a control 
panel) and in turn controls the various subsystems. The 
host computer 20 performs system level control func- 
5 tions. It accepts inputs from the operator via the operator 
interface 22 as well as system status changes (e.g., 
mode changes) and makes appropriate system chang- 
es. A system control bus (not shown) provides the inter- 
face from the host computer to the subsystems. A scan 
10 controller (not shown) provides real-time (acoustic vec- 
tor rate) control inputs to the various subsystems. The 
scan controller is programmed by the host computer 
with the vector sequences and synchronization options 
for acoustic frame acquisitions. Thus, the scan control- 
is ler controls the beam distribution and the beam density. 
The scan controller transmits the beam parameters de- 
fined by the host computer to the subsystems via a scan 
control bus (not shown). 

[0031 ] The conventional system has the capability to 

20 superimpose graphical symbols on any ultrasound im- 
age. The superimposition of graphics on the image 
frame is accomplished in the video processor 14, which 
receives the ultrasound image frame from the X-Y dis- 
play memory in the scan converter 12 and the graphics 

25 data from graphics/timeline display memory 18. The 
graphics data is processed and input into the graphics/ 
timeline display memory 1 8 by the host computer 20 or, 
alternatively, by a dedicated graphics processor which 
is synchronized with the other subsystems by the host 

30 computer. The host computer is programmed to monitor 
the position of a trackball manipulated by the system op- 
erator on the operator interface, acquire spectral Dop- 
pler imaging data from a sample volume determined by 
the trackball position, and superimpose a sample gate 

35 cursor on the displayed image frame at a location cor- 
responding to the trackball position. Similarly, the host 
computer is programmed to monitorthe state of a toggle 
switch on the operator interface and to control the size 
of the sample volume (and the sample gate cursor) as 

40 a function of the state of the toggle switch. 

[0032] In accordance with the preferred embodiment 
of the invention, the Doppler sample gate-position and 
size, as well as the beam steering angle, are automati- 
cally set based on actual vessel image data. This auto- 

45 matic adjustment can be performed either when the 
sample gate cursor is first activated or brought up in the 
image (before or upon spectral Doppler mode activa- 
tion); or when the sample gate cursor is being moved 
by the system operator. FIG. 2 represents an image 

so frame 32 produced by the system shown in FIG. 1 . The 
exemplary image shown in frame 32 includes a visual 
representation of a blood vessel 30 with conventional 
graphics superimposed thereon. The displayed graph- 
ics include: a Doppler beam cursor (beam centerline) 

55 26; a vessel slope cursor 28; and Doppler sample gate 
(sample volume) graphic, consisting of a sample gate 
top graphic 36 and a sample gate bottom graphic 34. In 
the latter type of sample gate graphic, the size of the 
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sample gate is represented by the distance separating 
the bottom graphic 34 and the top graphic 36 in FIG. 2. 
However, it will b e readily app reciated by persons skilled 
in the art that the sample gate graphic may have a dif- 
ferent geometry, in which case the size of the sample 
gate would correspond to a dimension of such graphic, 
e.g., a diameter if the graphic were a circle. The esti- 
mated value of the Doppler angle between the Doppler 
beam cursor 26 and the vessel slope cursor 28 on the 
vessel 30 in the image 32 is used to convert Doppler 
frequency shifts into velocity units according to Eq, (1). 
The Doppler angle value is usually displayed along with 
the graphic. 

[0033] In accordance with the preferred embodiment 
of the invention, the position and size of the sample vol- 
ume (and the sample gate cursor formed by graphics 34 
and 36) are set automatically by the ultrasound imaging 
system. In addition, the ultrasound imaging system can 
be programmed to select a steering angle which mini- 
mizes the Doppler angle, i.e., the angle between the 
beam steering direction and the vessel slope cursor. 
[0034] To set the sample volume or gate position au- 
tomatically, the system first searches for the best vessel 
segment within the target search region. The target 
search region depends on the current Doppler operating 
condition. For the case when the sample gate cursor is 
first activated or brought up in the image (before or upon 
activation of the spectral Doppler imaging mode), the 
target search region may be a predetermined region of 
interest centered at a preset sample gate position. For 
the case when the sample gate cursor is being moved, 
the image quadrant (or a smaller sector) that the sample 
gate cursor is moving towards can be used to define the 
search region. For example, FIG. 3 shows a rectangular 
sector 38 in the upper left-hand quadrant of the image 
frame, toward which the sample gate cursor is moving. 
The image frame 32 shows portions of three vessels 30, 
30' and 30". In accordance with one aspect of the pre- 
ferred embodiment, the search region 38 is searched to 
find the best vessel segment therein. In this example, 
the best vessel segment, over which the re-selected 
sample gate is positioned, is a segment of vessel 30'. 
The re-selected sample gate is represented on the dis- 
play screen by the sample gate bottom graphic 34' and 
the sample gate top graphic 36'. 
[0035] In general, a target search region may contain 
more than one vessel segment. The "best" vessel seg- 
ment can be chosen based on any reasonable morpho- 
logical characteristic, such as the vessel diameter (clos- 
est to a standard size of the vessel type associated with 
the user-selected application type), greatest vessel 
length or area, or most uniform diameter (since the user 
usually moves the probe to obtain the best view of the 
vessel of interest), or a combination of goodness meas- 
ures. Alter-natively, the best vessel segment can be de- 
fined as the one which is at the shortest distance from 
the preset position (for the case when the sample gate 
cursor is first activated) or current sample gate position 



(for the case when the sample gate cursor is being 
moved). 

[0036] To identify the best vessel segment within the 
search region, the algorithm shown in FIG. 4 is used. In 
5 a typical scanner, the current image frame can be read 
(step 40) either from the X-Y display memory in the scan 
converter or from the cine memory. To facilitate the ves- 
sel segment detection process, the image is first subject 
to some preprocessing (step 42), which may include any 
10 or all of the following: (1) image decimation to reduce 
the image to a smaller two-dimensional array of pixel 
elements; (2) def in ition of a search region within the dec- 
imated image; (3) application of a smoothing filter (e.g., 
median or mean filter) to reduce speckle noise; and (4) 
15 pixel intensity histogram equalization to enhance con- 
trast between background and vessel ("contrast stretch" 
in block 42 in FIG. 4). The image decimation factor (e. 
g., 2) is in general dependent on a predefined minimum 
vessel size and the image size (in terms of pixel count) 
for the current image depth setting. 
[0037] Still referring to FIG. 4, an adaptive threshold- 
ing method is used to binarize the search region (step 
46) during B-mode imaging. The objective is to segment 
the search region by marking each pixel as either "1 " if 
it is inside a cystic structure or "0" if it is not. Assuming 
that blood flow is hypoechoic, whereas soft tissues have 
relatively great pixel intensities, segmentation is 
achieved by comparing the intensity of each pixel to a 
threshold. In the simplest implementation, the threshold 
may be chosen based on some fraction of the global 
maximum intensity in the search region. For more robust 
performance, however, an adaptive scheme is also pro- 
posed wherein a threshold based on the local pixel in- 
tensity distribution is used to classify each pixel. 
[0038] While many variations of adaptive thresholding 
schemes are possible, a preferred binarization method 
for B-mode imaging is shown in FIG. 5. For each new 
pixel Py in the search region (step 56), we first compute 
the local mean (M) and standard deviation (a) of the 
neighboring pixels within an RxR kernel centered at pix- 
el Pjj (step 58). The kernel size R (e.g., 9) is chosen to 
be a little smaller than the expected minimum vessel 
size. 

[0039] To test if the kernel contains only tissue speck- 
le and is not inside a vessel, the standard deviation o is 
compared to a pre-defined threshold value $ (step 60). 
If g><|>, the pixels in the kernel are highly unlikely to be 
"blood flow"; however, if the kernel is partially inside a 
vessel, the condition of o>$ may still hold. Therefore, 
only pixels in the kernel with a value greater than the 
mean (M) are marked as "0" (step 62). If a<4>, an adap- 
tive threshold 7~ is computed (step 64) as follows: 

T=M-kXc (2) 

where k is a constant. In general, k is a small positive 
fraction such as 0.2. 
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[0040] If Pg <T (step 66 in FIG. 5), then P i} is marked 
as "0" (not blood flow) (step 68); otherwise it is marked 
as "1 " (blood flow) (step 70). The same marking step is 
repeated for each unmarked pixel in the search region 
(return to step 56). 

[0041] For best performance, both B-mode and color 
flow (velocity or power) image data should be used. If 
color mode is not turned on, a hidden color frame can 
be fired just for this purpose. The advantages of using 
a color image are that it usually represents a clearer im- 
age of vessels with flowing blood, and it is already seg- 
mented to highlight the vessel lumen with flow. That is, 
the colorized pixels can be directly treated as the Vs in 
a binarized image mask (step 44 in FIG. 4). The rest of 
the processing (object counting and characterization) is 
then the same as for B-mode data. Depending on which 
set of data provides more clear image segmentation 
(more clear images of vessel segments), the decisive 
vessel morphology or distance measurements can be 
based on either or both of the B-mode and color flow 
data. 

[0042] In general, the binarized image after adaptive 
thresholding may be very- "noisy," consisting of many 
small, isolated structures. These small structures (most- 
ly speckle noise) can be eliminated by using morpho- 
logical filters (step 48 in FIG. 4), which are nonlinear im- 
age transformation techniques taught in many digital im- 
age processing textbooks (see, e.g., William K. Pratt, 
Digital Image Processing, 2nd edition, Wiley, New York). 
Basic morphological operations can be implemented by 
direct pattern-matching ("hit or miss") transformations, 
or by using a more efficient pixel stacker and look-up 
table method. 

[0043] Erosion and dilation represent two basic mor- 
phological operations which, when used in series, can 
be quite effective in closing up the speckle noise struc- 
tures. Basically, each pass of an erosion filter strips off 
the outermost pixel layer of a continuous bright ("1 ") re- 
gion. This tends to close up the small extraneous bright 
structures, like speckle noise. The erosion operation will 
also erode the outermost layer of any blood flow region. 
To offset this undesirable effect, an opposite operation, 
called a dilation, can be applied after each pass of an 
erosion filter. The effect of a dilation filter is to add back 
a layer of pixels to existing bright objects. Speckle noise 
gaps which have been completely closed up (no longer 
exist) by erosion filtering will not be regenerated by the 
dilation filter. In practice, one pass of the erosion filter 
followed by a dilation filter can eliminate a majority of 
speckle noise structures. But if necessary, additional 
passes of erosion and dilation can be performed. 
[0044] The aim of step 60 in FIG. 4 is to count the total 
number of continuous objects in the morphologically fil- 
tered search region. A continuous object is defined as 
a set of object pixels ("1 M s) that are connected to each 
other by pixels within the set. In general, a continuous 
object may contain holes or rough boundary pixel layers 
with small extrusions. These irregularities can cause 



confusion during the object segmentation process. 
Therefore, in accordance with the preferred embodi- 
ment of the invention, it is assumed that the objects of 
interest (vessels) have relatively well-defined shapes 

5 without extremely large holes or extrusions. The objec- 
tive is to segment or group ail the object pixels ("1"s) 
into distinct objects with different object numbers. 
[0045] A computationally efficient approach that 
should work well in most situations is the following. Sup- 
to pose P i} denotes the binarized pixel value in row / and 
column j. Beginning from one corner of the search re- 
gion (e.g., P o0 ), the pixels in the image are checked and 
assigned an object number in a sequential manner as 
in a raster scan (i.e., row by row). If P,y = 0 (not blood 

15 flow), then it remains as zero. If P# = 1 (blood flow), a 
predefined kernel of neighboring pixels, which have al- 
ready been scanned, is reviewed. If at least one pixel in 
this neighborhood kernel is numbered as an object, then 
the current pixel is considered part of the same object 

20 and is assigned the same number. Otherwise, the cur- 
rent pixel is considered the beginning of a new object 
and is assigned a new object number. The neighbor- 
hood kernel consists of pixels in the previous column 
and/or previous row (which have already been 

25 scanned), i.e., pixels in row /, columns /-1 , y-2, y-3, .... /- 
n, and row /-1 , columns y-1 , y, /f 1 , J+2 t y+m where n 
and m are integer values. At a minimum, the pattern 
should include the three immediate neighbors P^,, 
P M ^ lf P Mj which have already been scanned. The 

30 checking of additional pixels in the previous row/col- 
umns is recommended to help distinguish between 
small extrusions and/or broken fragments of an object 
already numbered and the beginning of a new object. 
[0046] Having assigned numbers to all of the objects, 

35 the total number of objects can be counted easily. Also, 
the area of each object can be computed simply by 
counting the number of the pixels marked with a given 
object number. 

[0047] Still referring to FIG. 4, in step 52 objects that 

40 are too large or too small are screened out. Although in 
step 48 morphological filters were used to remove small 
speckle-like structures, some of them may still be 
present in the search region. These small structures will 
now have unique object numbers, but they can be elim- 

45 inated by checking whether their area is within a prede- 
fined range. If the area is smaller than some lower 
threshold (based on the mean speckle size of the scan- 
ner for the current imaging setup), it will be rejected. For 
example, a two-pixel object (Area = 2 pixels) represent- 

50 jng residual speckle noise will be rejected if the lower 
threshold is set at Area = 3 pixels. If the counted area 
of an object is too large to be the interior of a vessel, the 
object will also be rejected. Only the remaining objects 
are considered to be "true" blood flow. 

55 [0048] Following the rejection of objects having a size 
outside predetermined limits, the system automatically 
determines the "best" vessel segment in accordance 
with one of the goodness measures previously de- 
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scribed and then calculates the center point of that ves- 
sel segment. Although not shown in FIGS. 2 and 3, a 
curved blood vessel, when intersected by the scan 
plane, may appear in the image as a discrete object 
within the frame boundaries, rather than as an object 
which traverses the frame from one edge to another. In 
the either case, the center point can be determined in 
various ways. For example, in step 54 (see FIG. 4) the 
x and y coordinates of the vessel center can be comput- 
ed by simply taking the average values of the vessel 
boundaries along the horizontal and vertical axes, re- 
spectively. Other measures such as center-of-mass can 
also be used. For example, the "center of mass" of each 
remaining object could be computed using the following 
equations: 



2>, 



x n = 



- JsL 



N 



(3) 



2>. 



(4) 



[0049] In Eqs. (3) and (4), x 0 and y 0 are the coordi- 
nates of the center of mass of the object, N is the number 
of the pixels comprising the object, and X/and //are the 
row and column coordinates for pixel Pg. 
[0050] Having identified the center of the best vessel 
segment, the vessel diameter and orientation can be 
computed using any automated method. For example, 
the automated methods for measuring the Doppler an- 
gle taught in either U.S. Patent No. 5,690,116 or 
6,068,598 could be used. The method taught in the latter 
prior art reference will be described here with reference 
to FIGS. 6-8. 

[0051] Referring to FIG. 6, the data representing the 
center point of the "best" vessel segment is identified by 
the host computer as the center point 72 of a search 
region 76. Based on the B-mode and color flow image 
data found within the search area 76, a Doppler angle 
is computed. 

[0052] Referring to FIG. 8, after the center point of the 
search region has been identified (step 84), the host 
computer determines (step 86) whether the image 
frame includes color flow data at pixel addresses corre- 
sponding to the location of the center point. If the image 
frame includes color flow data corresponding to the 
center point, then the host computer searches out from 
the center point 72 along radial lines 74 which are an- 
gularly spaced S degrees apart over an entire 360° 
range, as shown in FIG. 6. The distance to search from 



the center is D cm. This edge search area is indicated 
by circle 76 in FIGS. 6 and 7. 

[0053] Along each radial line 74, the host computer 
searches from the center point 72 and stores the point 
5 as an edge point if it is the first of X points displaying B- 
mode intensity data instead of color flow velocity orpow- 
er data (step 88 in FIG. 8). Exemplary edge points are 
depicted as rectangles 78 in FIG. 7. If no such point is 
found before searching D cm or finding the edge of the 

10 color region of interest, then no edge point is marked 
along that radial line. Once each radial line is searched, 
all the edge points 78 in a certain segment of the edge 
point search area (e.g., Segment 1 indicated by dotted 
lines in FIG. 7) are grouped together (step 90) and fed 

15 to a line-fit algorithm which generates both a vessel 
slope estimate and a goodness-of-fit measurement 
(step 92). This is repeated for other segments (e.g., 
Segment 2 indicated by solid straight lines in FIG. 7) and 
in each case the slope of the vessel and the goodness 

20 of fit are recorded. The segments may overlap each oth- 
er by some number of degrees, as do Segments 1 and 
2 shown in FIG. 7. If a particular segment does not have 
some minimum number of edge points within it, then that 
segment is disregarded. 

25 [0054] In addition to the foregoing, the algorithm also 
determines B-mode edge points (step 94) by searching 
the B-mode intensity data from the center point in radial 
lines spaced S degrees apart over an entire 360° range. 
The distance to search from the center is D cm. Along 

30 each radial line each B-mode intensity value (corre- 
sponding to respective pixels) is replaced with the aver- 
age of itself and its two neighbors along the radius. The 
peak and minimum intensities along the averaged radial 
line as well as the largest difference (from one pixel to 

35 the next) are each recorded. If the difference between 
the peak and minimum intensities does not exceed 
some threshold, then no edge point is specified for this 
ray. If the difference between the peak and minimum in- 
tensities exceeds the threshold, then a search is started 

40 at some number of points from the center and stops 
when a point (the edge point) is found to exceed a dif- 
ference-only threshold, an intensity-only threshold or a 
combined difference and intensity threshold. For exam- 
ple, if the pixel location is 50% of the maximum intensity 

45 and 30% of the maximum difference, then it would pass 
the combined difference and intensity threshold. The in- 
tensity at the edge point is noted. If no such point is 
found before searching Don or finding the edge of the 
B-mode image, then no edge point is marked along that 

50 radial line. Once each radial line has been searched, 
some percent of the edge points are disregarded. The 
disregarded edge points are those associated with the 
lowest intensities. All of the remaining edge points in a 
certain segment of the edge point search area are 

55 grouped (step 96) and then fed to a line-fit algorithm 
which generates both a vessel slope estimate and a 
goodness-of-fit measurement (step 98). This is repeat- 
ed for other segments, and in each case the vessel 
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slope and the goodness of fit are recorded. The seg- 
ments may overlap each other by some number of de- 
grees. If a particular segment does not have some min- 
imum numberof edge points within it, then that segment 
is disregarded. 

[0055] If no B-rnode or color flow mode segment gen- 
erated enough edge points to get a vessel slope esti- 
mate, the distance D is increased and the algorithm is 
rerun. 

[0056] At this point in the algorithm, estimates of the 
vessel slope and their corresponding goodness-of-fit 
measurements are known for some number of seg- 
ments (for B mode and color flow mode). The segment 
having the best goodness of fit is identified (step 100) 
and its vessel slope is combined (averaged) with all the 
other vessel slope estimates that have a goodness of fit 
measurement not exceeding some difference relative to 
the best one (step 1 02). However, if color is active and 
the best color vessel slope exceeds some number of 
degrees (indicating a vessel somewhat vertical), then 
only color data is used in this vessel slope-combining 
algorithm. This is done because somewhat vertical ves- 
sel walls are difficult to detect in B-mode due to the lat- 
eral smearing of the image data. The host computer de- 
termines the vessel diameter (step 1 04 in FIG. 8) by cal- 
culating the distance between edge points along a line 
passing through the center point and orthogonal to the 
calculated vessel slope. 

[0057] After the vessel diameter and orientation have 
been computed, the sample gate size can be automat- 
ically set (step 1 06 in FIG. 8) to be some fraction of the 
vessel diameter. The optimal fraction can be pre-select- 
ed based on the practice of the individual clinical site. It 
does not have to be 1 00% or greater, since the typical 
Doppler exam is aimed at detecting arterial diseases 
and not for volume flow measurement. 
[0058] Alternatively, if power Doppler image of the 
vessel is available, the sample gate size adjustment can 
be based on power Doppler intensity levels within the 
vessel lumen. This technique makes use of the fact that 
the power Doppler intensity usually varies from zero 
near the vessel wail to maximum brightness towards the 
center of the vessel. For example, the sample gate 
boundaries can be placed at the -1 0 dB drop-off points 
from the maximum in the power Doppler image of the 
vessel lumen. 

[0059] Finally, having established the vessel orienta- 
tion based on image analysis, the best steering angle 
can be chosen automatically (step 1 08 in FIG. 8) to min- 
imize the Doppler angle. For example, if the imaging 
system is programmed with multiple sets of beam-form- 
ing time delay tables corresponding to multiple beam 
steering angles, the host computer need only determine 
which beam steering angle is closest to the vessel slope 
calculated based on image analysis, i.e., by identifying 
which beam steering angle minimizes the difference be- 
tween the beam steering angle and the vessel orienta- 
tion angle. 



[0060] Although the preferred embodiment has been 
disclosed in the context of an ultrasound imaging sys- 
tem having a host computer which processes graphics, 
a separate dedicated graphics processor can be em- 
5 ployed which is separate from the host computer. Fur- 
thermore, although FIG. 1 shows the architecture of an 
ultrasound imaging system having multiple processors 
and a host computer, all of the data processing and com- 
puting functions could be performed by a single compu- 
te ter having sufficient processing power. 

[0061] For the sake of good order, various aspects of 
the invention are set out in the following clauses:- 

1 . A method of locating a sample gate in a test sub- 
's ject, comprising the steps of: 

(a) detecting a predetermined state associated 
with a sample gate graphic superimposed on 
an image being displayed; 

20 (b) processing a frame of imaging data from 

which said image was derived to determine a 
predetermined point in a vessel segment in said 
test subject in response to detection of said pre- 
determined state; 

25 (c) moving a sample gate to a predetermined 

position relative to said predetermined point in 
response to completion of said processing 
step; and 

(d) superimposing a sample gate graphic at a 
30 new position in said image, said new position 

corresponding to said predetermined position 
of said sample gate. 

2. The method as recited in clause 1 , further com- 
35 prising the steps of: 

(e) processing said imaging data to determine 
a dimension of said vessel segment; and 

(f) adjusting a dimension of said sample gate 
40 to have a predetermined relationship to said di- 
mension of said vessel segment. 

3. The method as recited in clause 2, wherein said 
dimension represents an interna! diameter of said 

45 vessel segment. 

4. The method as recited in clause 1 , wherein said 
predetermined state comprises activation of said 
sample gate graphic. 

50 

5. The method as recited in clause 1 , wherein said 
predetermined state comprises movement of said 
sample gate graphic. 

55 6. The method as recited in clause 1 , wherein said 
step (b) comprises the steps of: 

defining a search region in said image; and 
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searching said search region for a continuous 
object representing a vessel segment having a 
predetermined morphological characteristic, 
wherein said sample gate graphics is placed in 
predetermined relationship to a center point of 5 
said continuous object representing said vessel 
segment having said predetermined morpho- 
logical characteristic. 

7. The method as recited in clause 6, wherein said 10 
searching step comprises the step of binarizing the 
pixel values within said search region. 

8. The method as recited in clause 7, wherein said 
binarizing step comprises the step of binarizing the is 
pixel values based on a pixel value threshold that 

is adaptive to local pixel value statistics. 

9. The method as recited in clause 7, wherein said 
binarizing step comprises the step of binarizing the 20 
pixel values based on the presence versus absence 

of color flow information in each pixel. 



1 5. The method as recited in clause 2, wherein said 
step (e) is based on power Doppler intensity levels 
acquired for pixels within said vessel segment. 

1 6. A method of sizing a sample gate in a test sub- 
ject, comprising the steps of: 

(a) detecting a predetermined state associated 
with a sample gate graphic being displayed; 

(b) processing a frame of imaging data from 
which said image was derived to determine a 
dimension of a vessel segment in said test sub- 
ject in response to detection of said predeter- 
mined state; 

(c) adjusting a dimension of said sample gate 
to have a predetermined relationship to said di- 
mension of said vessel segment in response to 
completion of said processing step; and 

(d) superimposing a sample gate graphic on 
said image, said sample gate graphic having a 
dimension corresponding to said adjusted di- 
mension of said sample gate. 



1 0. The method as recited in clause 7, wherein said 
searching step further comprises the step of per- 25 
forming morphological filtering to eliminate struc- 
tures smaller than a speckle size. 

1 1 . The method as recited in clause 1 0, wherein 
said searching step further comprises the steps of: 30 

counting the number of continuous objects in 
said search region; 

computing the area of each continuous object; 
and 35 
rejecting objects having an area lying outside a 
predetermined range. 

12. The method as recited in clause 1 , further com- 
prising the steps of: *o 

/(d) processing said imaging data to determine 
an angle of orientation of said vessel segment 
in said image; and 

(e) selecting a beam steering angle from a plu- 
rality of predefined beam steering angles which 
minimizes an angle between said angle of ori- 
entation of said vessel segment and said beam 
steering angle. 

50 

1 3. The method as recited in clause 1 , wherein step 
(b) comprises the step of finding a center of mass 
of said vessel segment. 

1 4. The method as recited in clause 1 , wherein step 55 
(b) comprises the step of calculating average val- 
ues of a boundary of said vessel segment along hor- 
izontal and vertical axes. 



17. The method as recited in clause 16, wherein 
said dimension of said continuous object repre- 
sents an internal diameter of said vessel segment. 

18. A system comprising: 

a display device comprising a multiplicity of pix- 
els; 

a memory for storing a frame of image pixel val- 
ues; and 

a computer programmed to perform the steps 
of: 

(a) controlling said display device to dis- 
play an image derived from said frame of 
image pixel values and a sample gate 
graphic superimposed on said image; 

(b) detecting a predetermined state asso- 
ciated with said sample gate graphic being 
displayed; 

(c) processing said image pixel values to 
determine a predetermined point in a ves- 
sel segment in said test subject in re- 
sponse to detection of said predetermined 
state; 

(d) moving a sample gate to a predeter- 
mined position relative to said predeter- 
mined point in response to completion of 
said processing step; and 

(e) superimposing a sample gate graphic 
at a new position in said image, said new 
position corresponding to said predeter- 
mined position of said sample gate. 

1 9. The system as recited in clause 1 8, wherein said 
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computer is further programmed to perform the fol- 
lowing steps: 

(f) processing said image pixel values to deter- 
mine a dimension of said vessel segment; and 5 

(g) adjusting a dimension of said sample gate 
graphic so that it has a predetermined relation- 
ship to said dimension of said vessel segment. 

20. The system as recited in clause 1 9, wherein said io 
dimension represents an internal diameter of said 
vessel segment. 

21 . The system as recited in clause 1 8, wherein said 
predetermined state comprises activation of said ^ 
sample gate graphic. 

22. The system as recited in clause 1 8, wherein said 
predetermined state comprises movement of said 
sample gate graphic. 2° 

23. The system as recited in clause 1 8, wherein said 
computer is further programmed to perform the fol- 
lowing steps: 

25 

defining a search region in said image; and 
searching said search region for a continuous 
object representing a vessel segment having a 
predetermined morphological characteristic, 
wherein said sample gate graphics is placed in 30 
predetermined relationship to a center point of 
said continuous object represent ng said vessel 
segment having said predetermined morpho- 
logical characteristic. 

35 

24. The system as recited in clause 23, wherein said 
searching step comprises the step of binarizing the 
pixel values within said search region. 

25. The system as recited in clause 24, wherein said 40 
binarizing step comprises the step of binarizing the 
pixel values based on a pixel value threshold that 

is adaptive to local pixel value statistics. 

26. The system as recited in clause 24, wherein said 
binarizing step comprises the step of binarizing the 
pixel values based on the presence versus absence 
of color flow information in each pixel. 

27. The system as recited in clause 24, wherein said so 
searching step further comprises the step of per- 
forming morphological filtering to eliminate struc- 
tures smaller than a speckle size. 

28. The system as recited in clause 27, wherein said 55 
searching step further comprises the steps of: 

counting the number of continuous objects in 



said search region; 

computing the area of each continuous object; 
and 

rejecting objects having an area lying outside a 
predetermined range. 

29. The system as recited in clause 1 8, wherein said 
computer is further programmed to perform the fol- 
lowing steps: 

(f) processing said imaging data to determine 
an angle of orientation of said vessel segment; 
and 

(g) selecting a beam steering angle from a plu- 
rality of predefined beam steering angles which 
minimizes an angle between said angle of ori- 
entation of said vessel segment and said beam 
steering angle. 

30. The system as recited in clause 18, wherein 
step (c) comprises the step of finding a center of 
mass of said vessel segment. 

31. The system as recited in clause 18, wherein 
step (b) comprises the step of calculating average 
values of a boundary of said vessel segment along 
horizontal and vertical axes. 

32. The system as recited in clause 1 9, wherein said 
step (f) is based on power Doppler intensity levels 
acquired for pixels within said vessel segment. 

33. The system as recited in clause 1 8, further com- 
prising: 

an ultrasound transducer array comprising a 
multiplicity of transducer elements; 
a transmit beam-former for pulsing selected 
transducer elements to transmit a series of ul- 
trasound transmit beams in a scan plane; 
a receive beam-former coupled to selected 
transducer elements of said transducer array 
for acquiring respective receive signals subse- 
quent to respective beam transmits; 
a signal processorforforming vectors of image 
parameter values from said receive signals; 
and 

a scan converter for converting said vectors in- 
to a frame of image pixel values and storing 
said frame of image pixel values in said mem- 
ory. 

34. A system comprising: 

a display device comprising a multiplicity of pix- 
els; 

a memory for storing a frame of image pixel val- 
ues; and 
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a computer programmed to perform the steps 
of: 

(a) controlling said display device to dis- 
play an image derived from said frame of 5 
image pixel values and a sample gate 
graphic superimposed on said image; 

(b) detecting a predetermined state asso- 
ciated with said sample gate graphic being 
displayed; 10 

(c) processing said image pixel values to 
determine a dimension of a vessel seg- 
ment in said test subject in response to de- 
tection of said predetermined state; 

(d) adjusting a dimension of said sample is 
gate to have a predetermined relationship 

to said dimension of said vessel segment 
in response to completion of said process- 
ing step; and 

(e) superimposing a sample gate graphic 20 
on said image, said sample gate graphic 
having a dimension corresponding to said 
adjusted dimension of said sample gate. 

35. The system as recited in clause 34, wherein said 25 
dimension represents an internal diameter of said 
vessel segment. 

36. A system comprising: 

30 

a display device; 

a memory for storing a frame of image pixel val- 
ues; 

means for controlling said display device to dis- 
play an image derived from said frame of image 35 
pixel values; 

means for controlling said display device to dis- 
play a sample gate graphic superimposed on 
said image; 

means for detecting a predetermined state as- *o 
sociated with said sampie gate graphic being 
displayed; 

means for processing said image pixel values 
to determine a predetermined point in a vessel 
segment in said test subject in response to de- *s 
tection of said predetermined state; 
means for moving a sample gate to a predeter- 
mined position relative to said predetermined 
point in response to completion of said process- 
ing step; and so 
means for controlling said display device to su- 
perimpose a sample gate graphic at a new po- 
sition in said image, said new position corre- 
sponding to said predetermined position of said 
sample gate. ss 

37. The system as recited in clause 36, further com- 
prising: 



means for processing said imaging data to de- 
termine a dimension of said vessel segment; 
and 

means for adjusting a dimension of said sample 
gate graphic to have a predetermined relation- 
ship to said dimension of said vessel segment. 



Claims 

1 . A method of locating a sample gate in a test subject, 
comprising the steps of: 

(a) detecting a predetermined state associated 
with a sample gate graphic superimposed on 
an image being displayed; 

(b) processing a frame of imaging data from 
which said image was derived to determine a 
predetermined point in a vessel segment in said 
test subject in response to detection of said pre- 
determined state; 

(c) moving a sample gate to a predetermined 
position relative to said predetermined point in 
response to completion of said processing 
step; and 

(d) superimposing a sample gate graphic at a 
new position in said image, said new position 
corresponding to said predetermined position 
of said sample gate. 

2. The method as recited in claim 1, further comprising 
the steps of: 

(e) processing said imaging data to determine 
a dimension of said vessel segment; and 

(f) adjusting a dimension of said sample gate 
to have a predetermined relationship to said di- 
mension of said vessel segment. 

3. A method of sizing a sample gate in a test subject, 
comprising the steps of: 

(a) detecting a predetermined state associated 
with a sample gate graphic being displayed; 

(b) processing a frame of imaging data from 
which said image was derived to determine a 
dimension of a vessel segment in said test sub- 
ject in response to detection of said predeter- 
mined state; 

(c) adjusting a dimension of said sample gate 
to have a predetermined relationship to said di- 
mension of said vessel segment in response to 
completion of said processing step; and 

(d) superimposing a sample gate graphic on 
said image, said sample gate graphic having a 
dimension corresponding to said adjusted di- 
mension of said sample gate. 
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The method as recited in claim 1 6 : wherein said di- 
mension of said continuous object represents an in- 
ternal diameter of said vessel segment. 

A system comprising: 5 

a display device (16) comprising a multiplicity 
of pixels; 

a memory (24) for storing a frame of image pixel 
values; and 10 
a computer (20) programmed to perform the 
steps of: 

(a) controlling said display device to dis- 
play an image derived from said frame of *5 
image pixel values and a sample gate 
graphic superimposed on said image; 

(b) detecting a predetermined state asso- 
ciated with said sample gate graphic being 
displayed; 20 

(c) processing said image pixel values to 
determine a predetermined point in a ves- 
sel segment in said test subject in re- 
sponse to detection of said predetermined 
state; 25 

(d) moving a sample gate to a predeter- 
mined position relative to said predeter- 
mined point in response to completion of 
said processing step; and 

(e) superimposing a sample gate graphic 30 
at a new position in said image, said new 
position corresponding to said predeter- 
mined position of said sample gate. 

The system as recited in claim 5, wherein said com- 35 
puter is further programmed to perform the follow- 
ing steps: 

(f) processing said image pixel values to deter- 
mine a dimension of said vessel segment; and 40 

(g) adjusting a dimension of said sample gate 
graphic so that it has a predetermined relation- 
ship to said dimension of said vessel segment. 

A system comprising: 45 

a display device (16) comprising a multiplicity 
of pixels; 

a memory (24) for storing a frame of image pixel 
values; and so 
a computer (20) programmed to perform the 
steps of: 

(a) controlling said display device to dis- 
play an image derived from said frame of 55 
image pixel values and a sample gate 
graphic superimposed on said image; 

(b) detecting a predetermined state asso- 



ciated with said sample gate graphic being 
displayed; 

(c) processing said image pixel values to 
determine a dimension of a vessel seg- 
ment in said test subject in response to de- 
tection of said predetermined state; 

(d) adjusting a dimension of said sample 
gate to have a predetermined relationship 
to said dimension of said vessel segment 
in response to completion of said process- 
ing step; and 

(e) superimposing a sample gate graphic 
on said image, said sample gate graphic 
having a dimension corresponding to said 
adjusted dimension of said sample gate. 

8. The system as recited in claim 6, wherein said di- 
mension represents an internal diameter of said 
vessel segment. 

9. A system comprising: 

a display device (16); 

a memory (24) for storing a frame of image pixel 
values; 

means (14, 20) for controlling said display de- 
vice to display an image derived from said 
frame of image pixel values; 
means (14, 18, 20) for controlling said display 
device to display a sample gate graphic super- 
imposed on said image; 
means (20) for detecting a predetermined state 
associated with said sample gate graphic being 
displayed; 

means (20) for processing said image pixel val- 
ues to determine a predetermined point in. a 
vessel segment in said test subject in response 
to detection of said predetermined state; 
means (20) for moving a sample gate to a pre- 
determined position relative to said predeter- 
mined point in response to completion of said 
processing step; and 

means (14, 18, 20) for controlling said display 
device to superimpose a sample gate graphic 
at a new position in said image, said new posi- 
tion corresponding to said predetermined posi- 
tion of said sample gate. 

10. The system as recited in claim 9, further comprising: 

means (20) for processing said imaging data to 
determine a dimension of said vessel segment; 
and 

means (20) for adjusting a dimension of said 
sample gate graphic to have a predetermined 
relationship to said dimension of said vessel 
segment. 
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B-MODE IMAGE 



BINARIZE SEARCH REGION 
BASED ON PRESENCE OR 
ABSENCE OF COLOR FLOW 
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BINARIZE SEARCH REGION 
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